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Abstract. Recently, attention has been drawn towards black 
carbon aerosols as a short-term climate warming mitigation 
candidate. However the global and regional impacts of the 
direct, indirect and semi-direct aerosol effects are highly un- 
certain, due to the complex nature of aerosol evolution and 
the way that mixed, aged aerosols interact with clouds and 
radiation. A detailed aerosol microphysical scheme, MA- 
TRIX, embedded within the GISS climate model is used in 
this study to present a quantitative assessment of the impact 
of microphysical processes involving black carbon, such as 
emission size distributions and optical properties on aerosol 
cloud activation and radiative effects. 

Our best estimate for net direct and indirect aerosol ra- 
diative flux change between 1750 and 2000 is —0.56 W/m 2 . 
However, the direct and indirect aerosol effects are quite sen- 
sitive to the black and organic carbon size distribution and 
consequential mixing state. The net radiative flux change 
can vary between —0.32 to —0.75 W/m 2 depending on these 
carbonaceous particle properties at emission. Taking into 
account internally mixed black carbon particles let us sim- 
ulate correct aerosol absorption. Absorption of black car- 
bon aerosols is amplified by sulfate and nitrate coatings and, 
even more strongly, by organic coatings. Black carbon mit- 
igation scenarios generally showed reduced radiative flux- 
eswhen sources with a large proportion of black carbon, 
such as diesel, are reduced; however reducing sources with a 
larger organic carbon component as well, such as bio-fuels, 
does not necessarily lead to a reduction in positive radiative 
flux. 
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1 Introduction 

Anthropogenic and natural aerosols impact the Earth’s radia- 
tion balance and thus exert a forcing on global climate. Black 
carbon (BC) has a positive radiative forcing, whileother 
aerosol species mainly lead to atmospheric cooling. There- 
fore cleaning up most aerosol emissions, which is much 
needed due to their severe health impacts, could strengthen 
the global warming trend. However, mitigation of BC has 
the potential to be beneficial in mitigating both climate 
warming and air pollution. Emission reductions that target 
light-absorbing aerosol might also reduce warming quickly 
(Hansen et al., 2000; Jacobson, 2002; Bond and Sun, 2005). 
However, BC’s indirect (cloud microphysics) and semi-direct 
contributions to net climate forcing remain an outstanding 
uncertainty. For example, the model study of Penner et 
al.. (2003) suggested that “smoke” did not produce net warm- 
ing, considering both direct and cloud effects. In addition, 
model experimental results are very sensitive to treatment of 
aerosol microphysics. 

The aerosol direct effect (ADE) is caused by absorption 
and scattering of solar radiation by liquid and solid aerosol 
particles in the atmosphere. Strongly absorbing, and there- 
fore enhancing climate warming, are black carbon (BC) par- 
ticles. The mass absorption strength and scattering ability of 
a BC particle strongly depends on the particle effective size 
and mixing state. Both of those quantities depend on the size 
and chemical composition of black carbon as it is released 
into the atmosphere as well as its microphysical and chemi- 
cal evolution during transport through the atmosphere. Most 
other aerosol species scatter radiation back to space, such 
as sulfate, nitrate, organic carbon (although some organics 
can be slightly absorbing), sea salt, aerosol water and to a 
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certain extent mineral dust, and therefore these counterbal- 
ance climate warming, by cooling the atmosphere. However, 
if BC is mixed with these scattering components, the result 
is enhanced absorption. Laboratory studies by Khalizov et 
al. (2009) showed that for soot particles with an initial mo- 
bility diameter of 320 nm and a 40% H 2 SO 4 mass coating 
fraction, absorption and scattering are increased by 1 .4 - and 
13 fold at 80% RH, respectively. Also, the single scattering 
albedo of soot aerosol increases from 0.1 to 0.5 after coat- 
ing and humidification. Jacobson [2000] performed global 
simulations with core- shell treatment and found BC forcings 
50% higher and 40% lower than forcings obtained with ex- 
ternally mixed and well-internally mixed treatments, respec- 
tively. Therefore, correct treatment of the aerosol size and 
mixing state is crucial in order to calculate their net radiative 
forcing. 

Attempts to quantify the aerosol direct radiative forcing, 
the difference in top of the atmosphere forcing between 
present day and pre-industrial conditions, by the AeroCom 
initiative gave estimates between —0.41 and +0.04 W/m 2 
[Schulz et al. 2006] and by the IPCC AR4 report a range 
from -0.9 to —0.1 W/m 2 [IPCC, 1007]. The wide range of 
uncertainty is caused on the one hand by large discrepancies 
among the model systems, ranging from emission strength, 
transport, aerosol transformation, removal and optical prop- 
erties. However, the models are also poorly constrained by 
limited availability of aerosol measurements. Furthermore, 
Myhre (2009) explains large deviations between models and 
satellite retrieved forcing estimates by the lack of consider- 
ing correct pre-industrial aerosol distributions in the satellite 
data derived forcing calculations and in failure to sample the 
model like the retrieval. 

The aerosol indirect effect (AIE) is caused by an increase 
in the number of hygroscopic aerosols, with a correspond- 
ing increase in cloud droplet number concentrations, reduc- 
tion in cloud droplet size, increase in cloud albedo (first in- 
direct effect) and suppression of precipitation and increased 
cloud life time (second indirect effect), thereby cooling the 
planet (IPCC, 2007). The role of BC in determining cloud 
droplet number concentrations (CDNC) is unknown [Adams 
and Seinfeld, 2003], as it depends upon the size distribu- 
tion and mixing state of BC with soluble species such as 
sulfate, nitrate and organic carbon. Furthermore, insoluble 
black carbon or dust aerosols can act as ice nuclei and there- 
fore alter cirrus and mixed-phase clouds; however these ef- 
fects are poorly known (e.g. Lohmann et al., 2008). Here 
again the microphysical evolution of the aerosol population 
in terms of size and mixing state is crucial in order to iden- 
tify the number of particles are suitable for cloud activation. 
The IPCC AR4 report estimates the aerosol indirect radiative 
forcing (AIF) to lie between —1.8 and —0.3 W/m 2 . The Ae- 
roCom initiative, linking models to satellite-based estimates 
gives a range of — 0.7±0.5 W/m 2 for the AIF[Quaas et al. 
2009]. The models on which these estimates were based 
generally did not include sophisticated treatment of aerosol 


mixing state, and the contribution of BC mixing to the AIF 
uncertainty was not evaluated. 

In this study we are particularly interested in the micro- 
physical evolution of black carbon. The size and mixing state 
of carbonaceous particles at emission must be assumed for 
this investigation. Once emitted, BC can grow by coagula- 
tion and condensation. Although freshly emitted pure BC is 
hydrophobic, inorganic and organic coatings will attract wa- 
ter and convert particles to be hygroscopic. Materials likely 
to condense on a BC particle, as on any other particle, are 
sulfate and nitrate precursors as well as secondary organics 
(Riemer et al. 2009). In addition to surface condensation 
processes, coagulation is a very effective mixing pathway. In 
the atmosphere every aerosol shape is unique, complex and 
not necessarily spherical. Transmission electron microscopy 
(TEM) images support the theory that BC particles become 
coated once emitted. Although BC may be internally mixed 
with other components, BC cannot be ‘well-mixed’ (diluted) 
in the particle, since soot, which contains BC, is irregularly 
shaped and mostly solid. Thus most likely BC is distinct and 
not well mixed within particles. 

As we explore the role of BC in aerosol forcing we will 
study the following questions: 

How important is the size distribution of carbonaceous 
particles at emission time for the overall aerosol simulation 
(Sect. 3)? How well can we constrain our model with obser- 
vational data (Sect. 4)? And how well can we assess BC re- 
duction strategies given the previously discussed microphys- 
ical modeling uncertainties (Sect. 5)? 

This paper is linked to a study by Menon et al. (2010) 
(hereafter referred to as SM10) where the same model ver- 
sion and configuration as in this study is used, but the analy- 
sis is more detailed with respect to the AIE. SM10 includes 
evaluation and comparison of the performance of anew AIE 
scheme coupled to aerosol microphysics or to the previous 
mass based GISS aerosol model (Koch et al., 2006; Menon 
and Rotstayn, 2006). 

2 Model description 

The Goddard Institute for Space Studies (GISS) General Cir- 
culation Model (GCM) climate modelE (Hansen et al., 2005; 
Schmidt et al., 2006) coupled to the aerosol microphysics 
and chemistry model MATRIX (Multiconfiguration Aerosol 
TRacker of mixing state) (Bauer et al., 2008), hereafter 
BA08, is used in this study. MATRIX is designed to sup- 
port model calculations of the direct and indirect effect and 
permit detailed treatment of aerosol mixing state, size and 
aerosol-cloud activation. 

MATRIX is based on the quadrature methods of moments 
(QMOM). MATRIX is under development to use the full ad- 
vantage of the QMOM scheme (McGraw 1997). In its cur- 
rent version MATRIX includes two moments, number and 
mass, and one quadrature point. Carrying only two moments 
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Table 1. Populations, constituents, and possible particle mixing state considerations in the radiation calculations: homogeneously internally 
mixed (HI), externally mixed (EM). 


population description 

symbol 

constituents 

mixing state options for the coupling 
between aerosols and radiation 

sulfate Aitken mode 

AKK 

SOJ 2 , noj,h 2 o 

EM/HI 

sulfate accum. mode 

ACC 

SO 4 2 , noj,h 2 o 

EM/HI 

dust accum. mode (<%5 inorg.) 

DD1 

dust, SOJ 2 , NO^,H 2 0 

EM 

dust accum. mode (>%5 inorg.) 

DS1 

dust.SO^ 2 , N07, H 2 0 

EM/HI 

dust coarse mode (<%5 inorg.) 

DD2 

dust.SO^ 2 , N07, H 2 0 

EM 

dust coarse mode (>%5 inorg.) 

DS2 

dust.SO^ 2 , N07, H 2 0 

EM/HI 

sea salt accum. mode 

SSA 

sea salt, SOJ 2 , N07, H 2 0 

EM/HI 

sea salt coarse mode 

SSC 

sea salt, SOJ 2 , N07, H 2 0 

EM/HI 

organic carbon (OC) 

OCC 

OC, SO^ 2 , NOJ, H 2 0 

EM/HI 

Black carbon (BC) (<5% inorg.) 

BC1 

bc,soj 2 ,no~,h 2 o 

EM/HI 

BC (5-20% inorg.) 

BC2 

BC, SO” 2 , NOJ, H 2 0 

EM/HI 

BC (>20% inorg.) 

BC3 

BC,SOJ 2 ,NO“,H 2 0 

HI 

BC-OC 

BOC 

BC, OC, S0” 2 ,N0“,H 2 0 

HI 

BC-sulfate (formed by coagulation) 

BCS 

BC,S04 2 ,N0“,H 2 0 

HI 

BC-mineral dust 

DBC 

dust, BC,S0^ 2 , NO^ , H 2 0 

HI 

mixed 

MXX 

dust, sea salt, BC, OC, SO^ 2 , NO 3 , H 2 0 

HI 


requires additional information about the shapes of the indi- 
vidual aerosol size distributions. We assume a lognormal dis- 
tribution with constant width when calculating the initial size 
distributions, the conversion between aerosol mass and num- 
ber concentration, emission distributions, coagulation rates 
and aerosol optical properties. 

For each aerosol population, defined by mixing state and 
size distribution, the tracked species are number concentra- 
tion, and mass concentration of sulfate, nitrate, ammonium, 
aerosol water, black carbon, organic carbon, mineral dust, 
and sea salt. Here we use the aerosol population setup called 
“mechanism 1” (BA08), given in Table 1. MATRIX dy- 
namics includes nucleation, new particle formation, particle 
emissions, gas-particle mass transfer, aerosol phase chem- 
istry, condensational growth, coagulation, and cloud activa- 
tion. 

In this paper the Napari et al. 2002 nucleation scheme 
is applied, based on classical nucleation theory involving 
H 2 SO 4 , NH 3 and H 2 O. The conversion from nucleation rate 
to new particles formation follows BA08. Coagulation can 
take place among all aerosol populations. Condensational 
growth is kinetically calculated for H 2 SO 4 deposition on all 
particle surfaces, and equilibrium is assumed between H 2 O, 
NH 3 and HNO 3 . The thermodynamic EQSAM model (Met- 
zger et al. 2006) is used to partition the total [NH 3 ]+[NH+ 4 ] 
and total [HN03]+[N0-3] between the gas and condensed 
phase, and to obtain the liquid aerosol water concentra- 
tion.Water uptake by sea salt aerosolis calculated after Lewis 
(2006). 


New additions to BA08 are linking the aerosol scheme 
to a recent aerosol indirect effect scheme (see Sect. 2.1)and 
the coupling of mixed aerosol populations to the radiation 
scheme as described in Sect. 2.2. 

2.1 Aerosol cloud coupling 

To simulate the indirect effect, we follow a similar treatment 
as described in Menon et al. (2010) that includes several 
changes to the treatment of cloud drop and ice crystal nu- 
cleation following the scheme from Morrison and Gettleman 
(2008). For cloud droplets, we use a prognostic equation to 
calculate CDNC, based on Morrison and Gettleman (2008) 
given as 

d(CDNC) _ 
dt 

Source— Loss (Autoconversion, Contact, Immersion) 

The source term includes newly nucleated CDNC and loss 
terms account for droplet loss through the process of au- 
toconversion (Qaut), contact nucleation and via immersion 
freezing. The source term is obtained from MATRIX us- 
ing the scheme of Abdul-Razak and Ghan (2000) that is 
based on the Kohler theory for multiple external lognormal 
modes that are composed of internally mixed soluble and in- 
soluble material. For ice crystal concentrations we include 
both heterogeneous freezing via immersion freezing based 
on Bigg (1953) and nucleation by deposition/condensation 
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freezing from Cooper (1986) following the treatment de- 
scribed in more detail by Morrison and Gettleman (2008). 
Here, aerosols do not directly affect ice crystal nucleation, 
except through their interaction with CDNC. The aerosol in- 
direct effect is obtained from the change in net (shortwave + 
longwave) cloud radiative forcing and may also be estimated 
from changes in net TOA radiation. 

2.2 Aerosol radiation coupling 

Previously the GISS radiation scheme treated only exter- 
nally mixed aerosol populations, suitable for our mass based 
aerosol scheme (Koch et al. 2006). Here we describe the 
new coupling scheme for internally mixed aerosol popula- 
tions. The GISS model radiation scheme (Hansen et al., 
1983) includes explicit multiple scattering calculations for 
solar radiation (shortwave - SW) and explicit integrations 
over both the SW and thermal (longwave - LW) spectral re- 
gions. Gaseous absorbers of LW radiation are H 2 O, CO 2 , O 3 , 
O 2 , and NO 2 . Size dependent scattering properties of clouds 
and aerosols are computed from Mie scattering, ray tracing, 
and T-matrix theory (Mishchenko et al. 1996) to include 
non- spherical cirrus and dust particles. The k-distribution 
approach (Lacis and Oinas 1991) utilizes 15 non-contiguous 
spectral intervals to model overlapping cloud aerosol and 
gaseous absorption. 

The MATRIX module calculates the optical properties sin- 
gle scattering albedo, asymmetry factor and extinction for the 
16 aerosol populations for six wavelength bands in the SW 
and 33 bands in LW and passes those to the GISS modelE ra- 
diation code. Each population can include multiple chemical 
species. If a particle contains only one chemical species, for 
example a freshly emitted pure organic carbon particle, then 
the optical properties are those of an externally mixed parti- 
cle. In those cases MATRIX calculates the refractive index 
for each population and then uses Mie-code pre-calculated 
lookup tables in order to assign the corresponding optical 
parameters. A homogeneous internal mixture is assumed 
to be a well-mixed particle, most likely containing soluble 
species including aerosol water. MATRIX calculates the op- 
tical properties for those particles by using the volume mix- 
ing approach. The following refractive indices at 550nm are 
used in this study: Sulfate ( 1.528- l.e-7/), nitrate (1.528-l.e- 
7/), OC (1.527-0.014/), BC (1.85-0.71/), sea salt (1.45-0./), 
dust (1.564-0.002/) and water (1.334-3.91e-8/). 

2.3 Model configuration 

Anthropogenic and natural emissions for present day and 
pre-industrial conditions are taken from the AeroCom project 
(Dentener et al., 2006 and http://nansen.ipsl.juissieu.fr/ 
AEROCOM/). We use fluxes for “natural” emissions of 
dust, sea salt and dimethyl sulfide (DMS), and organic car- 
bon (OC) assuming secondary organic aerosol as a 15% yield 
from terpene emissions, and in addition, anthropogenic emis- 
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sions from biomass burning and fossil and bio fuel burning 
of SO 2 , OC and BC. The inventory provides data for the year 
2000 (present-day conditions) and for the year 1750 (pre- 
industrial conditions). 

We will discuss direct and indirect aerosol effects by look- 
ing at aerosol direct radiative impacts (ADE) and cloud ra- 
diative impacts (AIE) at present-day conditions. ADE is 
calculated by taking the difference between radiative trans- 
fer calculations at the top of the atmosphere with and with- 
out aerosols. AIE is calculated from changes to the net 
cloud forcing obtained from differences between total and 
clear skies for each call to the radiation excluding aerosols. 
Aerosol direct (ADFC) and indirect radiative flux changes 
(AIFC) arecalculated as the difference in radiative fluxes be- 
tween pre-industrial and present day conditions. AIFC and 
AIE include indirect and semi-direct effects. Semi-direct ef- 
fects are cloud changes caused by changed vertical heating 
profiles in the atmosphere caused by aerosol forcing changes. 

All radiative flux calculations, ADE, AIE, ADFC and 
AIFC, are calculated in the aerosol-cloud coupled model, and 
therefore internal climate feedbacks impact all radiative flux 
calculations. Note that we do not use the terminology radia- 
tive forcing as this would require an offline forcing calcula- 
tion that does not contain climate feedbacks. 

Model simulations for present day and pre-industrial con- 
ditions only differ by emission levels. The model runs 
are performed fully interactive therefore aerosol and cloud 
changes will impact the general circulation. Historic sea 
surface temperatures are prescribed as boundary conditions. 
The model is not coupled to the gas-phase ozone chemistry 
scheme. Off-line fields of oxidants and nitric acid concen- 
trations are provided as input fields from a previous coupled 
chemistry - aerosol simulation. Off-line ozone fields are pro- 
vided for forcing calculation in the model radiation scheme. 
All offline chemistry fields and GHG concentrations are set 
to present day concentrations. Present day atmospheric con- 
ditions, including present day sea surface temperature, are 
also applied for the pre-industrial emission runs. The model 
is employed on a horizontal resolution of 4°x5° latitude by 
longitude and 23 vertical layers. The model uses a 30-min 
time step for all physics calculations. Every model simula- 
tion is integrated for 5.5 years, and if not otherwise noted, 
five-year mean conditions are discussed in this paper. 

3 Effects of chosen particle size for emission 

In order to calculate emission fluxes in microphysical aerosol 
models, the mass, number concentration and the mixing state 
of those emissions need to be known. Information about size 
distributions is very important for particulate emission fluxes 
such as dust, sea salt and carbonaceous emissions. Usually 
some size information is provided for dust and sea salt emis- 
sions, or interactively calculated in models, but not for car- 
bonaceous aerosols. Current emission inventories, such as 
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Bond et al. [2004] or Cooke et al. [1999] provide only mass 
emission information, so that each modeler has to choose 
what sizes to assign for carbonaceous emission fluxes. Tex- 
tor et al [2006] summarized the sizes of the emitted particles 
as used by 16 different models that participated in the Ae- 
roCom study. Mass median diameter ranging from 0.02 to 
0.85 pm were used for BC and OC emissions. Some Aero- 
Com models used different sizes for fossil and biomass burn- 
ing sources, and some models emitted BC and OC emissions 
into several size bins. These choices differed greatly among 
the models; furthermore there is very little information avail- 
able about actual emission sizes. Bond et al (2006, Table 
3.) collected particle size distribution observations at com- 
bustion sources from the literature and reported mass me- 
dian diameters of 0.038-0.32 pm for diesel vehicles, 0.02- 
1.5 pm for gasoline vehicles, 0 . 1 - 1 . 3 pm for small solid fuel 
combustors (excluding coal briquettes), such as wood fire- 
places or cooking stoves, and 0.05-0.78 pm for large station- 
ary sources such as industrial boilers. However these sizes 
obtained close to emitters may not be appropriate for global 
model gridbox-scale initial particle sizes. 

In order to understand the importance of emission size dis- 
tribution of carbonaceous particles, we performed a series 
of sensitivity studies where only the sizes of carbonaceous 
emissions are varied. The total mass of the emissions is 
identical between the experiments, but the emitted number 
concentration varies among the experiments according to the 
particle size assumptions. Note that those sizes are just the 
emission sizes, and that aged particle sizes are determined 
by the subsequent microphysical processes. In our simu- 
lations, carbonaceous aerosols grow by roughly a factor of 
three, through coagulation, condensation and water uptake. 
We assume that fossil and bio fuel BC and OC emissions 
enter the atmosphere as external mixtures, entering popula- 
tion BC1 and OCC, and biomass burning emissions are inter- 
nally mixed, entering population BOC in our model. How- 
ever the model is rather insensitive towards this last assump- 
tion, as BC and OC from biomass burning sources coagulate 
very quickly even if emitted as external mixtures (not shown 
here). 

The geometric mean emission particle diameters chosen 
for the single experiments are reported in Table 2. Three 
sensitivity experiments are performed around the base case 
simulation (BASE), one experiment with smaller emission 
sizes (SI) and two experiments with larger emission sizes 
(S3, S4).The set of experiments covers approximately the 
range of observed emission sizes. 

The design of the sensitivity experiments is rather simple, 
but analyzing the results is difficult as nearly everything in 
the aerosol simulation is affected. First we will analyze the 
differences between the base case BASE and S 1 , the exper- 
iment with the smallest BC/OC emission sizes. The largest 
differences occur in the coagulation calculation. Smaller car- 
bonaceous particles in experiment S 1 impact the mixing state 
and size distribution in all of MATRIX’S populations. Size is 
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Table 2. Emission sizes, geometric mean mass diameters in [/zm] 
of carbonaceous aerosols. 



OC [OCC] 
fossil & bio- 
fuel 

BC 

[BC1] 

fossil & bio- 
fuel 

BC-OC 

[BOC] 

biomass 

burning 

SI 

0.01 

0.01 

0.025 

BASE 

0.05 

0.05 

0.10 

S3 

0.1 

0.1 

0.25 

S4 

0.5 

0.5 

1 


the most important particle property affecting coagulation. 
The differences in sulfate, OC and BC mixing state between 
BASE and SI (SI-BASE) are presented in Fig. 1 for present 
day conditions. Smaller initial OC/BC sizes lead to more 
coagulation between sulfate and BC. The sulfate coagulation 
rate to form BCS (a BC-sulfate mixture) increases from 1087 
(BASE) to 8655 (SI) Gg/a, leaving less externally mixed 
sulfate in the system (sulfate ACC and AKK loads decrease 
from 188 (BASE) to 41 Gg (SI)). Table 3 in BA08 may be 
consulted to put those budget numbers into perspective. Net 
sulfate condensation rates are rather similar between the ex- 
periments, but the particles on which surfaces condensation 
happens changes. In BASE most H 2 SO 4 condenses on ACC, 
OCC and BOC particles and aqueous phase H 2 SO 4 conden- 
sations within clouds favors ACC and BOC. In SI most con- 
densation forms sulfates on BOC surfaces. Again leading to 
more internally mixed particles, due to less condensation of 
H 2 SO 4 on pure sulfate (ACC) particles. Changes in H2S04 
condensation rates also affect sulfate primary particle pro- 
duction, and decreases sulfate nucleation rates (1800 Gg/a 
(BASE), 1168 Gg/a (SI)), hence changing the size distribu- 
tion. BC (Fig. 1 row 2) itself favors coagulation with OC for 
SI, this however increases the sulfate ratio in the BCS popu- 
lation. So overall sulfate predominantly mixes with BC and 
BC mixes more with OC, leaving fewer externally mixed OC 
(OCC) in the atmosphere. Furthermore, the smaller overall 
OC/BC size leads to more coagulation of OC and BC with 
coarse particles. Figure 2 shows the zonal mean BC mass 
mixing ratio, the sum over BC from all BC containing popu- 
lations, of the BASE case and the differences between BASE 
and the size experiments SI, S3 and S4. The total BC mass 
increases from the BASE case 1195 Gg to 1287 Gg in SI, 
with all of the increase happening in the Northern Hemi- 
sphere. This is caused by the extended lifetime of smaller 
pollution BC particles and the fact that smaller particles can 
be transported over larger distances. 

Figure 3 presents the differences of experiments SI and 
S3 with BASE. Comparing BASE and SI, AOT increases 
over the oceans and decreases in biomass burning areas. The 
aerosol direct effect shows a similar pattern but of opposite 
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Fig. 1. Differences in present day total column aerosol mass concentrations [gg/m 2 ] per mixing state between experiment SI and BASE 
(SI-BASE). The first row shows externally mixed sulfate (population ACC and AKK), sulfate predominantly mixed with OC (OCC, BOC), 
BC (BC1, BC2, BC3, BCS), and sulfate mixed with coarse aerosols (DS1, DS2, SSC, MXX). Rows two and three show the corresponding 
changes in mixing state for BC and OC mixtures. Graphs show only the sulfate, or BC, OC respectively, mass present in those mixtures. 
Some aerosol concentrations are divided by a factor of ten to match the color bar as indicated in the individual title of the map. 


Table 3. Global annual average values for the net aerosol direct effect (ADE), for sulfate and nitrate (ACC), OC (OCC), BC, sulfate and 
nitrate (BC-SU, population BC1,BC2,BC3 BCS), BC, sulfate, nitrate and OC (BC-OC, population BOC); AIFC is indirect and semi-direct 
radiative flux change; ADFC aerosol direct radiatove flux change and NRtoa is top of the atmosphere net radiative forcing. 



ADEtoa 

ACC 

PD 

ADEtoa 

OCC 

PD 

ADEtoa 

BC-SU 

PD 

ADEtoa 

BC-OC 

PD 

ADEtoa 

total 

PD 

AIFC 

PD-PI 

ADFC 

PD-PI 

NRtoa 

PD-PI 

BASE 

-0.05 

-0.16 

0.05 

0.1 -1.78 

Emission size experiments 

-0.45 

-0.11 

-0.56 

SI 

-0.01 

-0.01 

0.01 

0.36 

-2.68 

-0.65 

-0.09 

-0.74 

S3 

-0.15 

-0.34 

0.09 

0.09 

-1.71 

-0.35 

-0.25 

-0.60 

S4 

-0.37 

-0.05 

0.03 

0.08 

-1.18 

-0.13 

-0.19 

-0.32 


sign. ADE increases over land, due to internal mixing of 
BC and therefore enhancing overall BC absorption, which 
leads to an increase in ADE due to BC particles. However 
the dominating feature here is the decreased ADE over the 
oceans, which is caused by a strong contribution from the 
mixed population, MXX. The smaller particles in SI lead 
to overall more mixing. If all 7 species are mixed together, 
or none of the other 15 populations can be used, MXX is 
populated or in other words MXX serves as population for 
all mixing states that can’t fit in any of the other 15 popu- 
lations. However MXX is always representing a coarse mode 
as usually some sea salt or dust is present in MXX. 


Changes in AIE depend on particle chemistry and num- 
ber concentrations, and therefore strongly on the size distri- 
bution, leading to different changes than seen in the mass 
concentrations as shown in Fig.l. Regionally the negative 
AIE (Fig. 3) is caused by the increased cloud droplet num- 
ber concentration (CDNC) in those regions. Most importan- 
tis the shift in size and mixing state. Particles are smaller, 
but still large enough to activate clouds, and more internally 
mixed in S 1 , leading to a larger number of cloud activating 
particles in polluted sulfate rich regions, but to regional de- 
creases in CDNCs in biomass burning, OC/BC rich regions. 
In the BASE case CDNCs in biomass burning regions are 
dominated by externally mixed OC, e.g. OCC, whereas in 
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Fig. 2. Zonal mean BC mass mixing [pg/kg] ratios for the BASE 
case and the differences between BASE (BA) and the size experi- 
ments, SI, S3, and S4. Global mean BC mass concentrations per 
experiment are given in the figure titles. 


SI more of OC is mixed with BOC leading to overall smaller 
number concentrations in biomass burning areas. Globally 
cloud radiative fluxes increase over the oceans leading to a 
difference of -0.83 W/m 2 in AIE. 

The second set of experiments, S3 and S4, tests the sen- 
sitivity of our model at larger BC/OC emission sizes. These 
experiments show many similarities towards each other, so 
that it is sufficient here to discuss the results of case S3. 
Larger initial OC/BC particles lead to less internal mixing 
of the aerosol populations, leaving more externally mixed 
sulfate and OC particles. Externally mixed sulfate increases 
from 188 Gg/a (BASE) to 575 Gg/a (S4) and more sulfate 
condenses on coarse particles. Also OC and BC are less 
mixed with each other, but BC in the biomass burning re- 
gions is more mixed with sulfate, leading again to an increase 
of sulfate in the BCS population. These changes in mixing 
state are generally the reverse of case S 1 ; however this does 
not translate into a reverse response for the aerosol optics 
(Fig. 3). AOT globally decreases due to less extinction by 
the slightly larger and more externally mixed particles. The 
direct aerosol effect (ADE) relative to BASE is slightly posi- 
tive, especially in dust regions. This is caused by less coagu- 
lation with the overall larger OC/BC/SU particles with dust, 
hence increasing dust lifetime in S4. The BC load (Fig. 2) is 
strongly reduced in S3 to 1041 Gg and 822 Gg in S4. 

CDNC particle concentrations are reduced by -15% (S3) 
and -20% (S4) leading to a weakened indirect effect, the 
aerosol indirect effect is approximated from changes to net 
radiation at the TOA and therefore also includes semi-direct 
effects, by 3% (S3) and 4% (S4); however there was also 
an increase in liquid water path (LWP) (0.5% (S3) and 0.8% 
(S4)) and increased cloud cover (1% (S3) and 0.3% (S4)). 


These sensitivity studies showed some surprising sensitiv- 
itiesbetween initial carbonaceous particle sizes and CDNC 
concentrations. The following processes cause CDNC 
changes: 

- Organic carbon is a hygroscopic, potentially cloud ac- 
tivating particle, e.g. assuming smaller/larger particle 
sizes at emission time will increase the emitted num- 
ber concentration and eventually lead to higher /lower 
CDNC. 

- Different ambient particle sizes, caused by the differ- 
ent assumption of initial carbonaceous aerosol sizes, 
impact as well for example sulfate containing parti- 
cles due to the changes in condensation and coagulation 
rates among the aerosol populations and impact primary 
aerosol formation due to the competition for precursor 
species between the processes of primary aerosol for- 
mation and aerosol growth or aging processes. 

After discussing the complex interactions a simple sum- 
mary (Fig. 4) can be given for the impact of carbonaceous 
aerosol emission sizes on radiative fluxes: Smaller BC/OC 
emission sizes lead to larger, more negative direct and in- 
direct aerosol effectsand larger BC/OC emission sizes to 
smaller ADE and AIE effects. 


4 Observational constraints 

The experiments performed in this study show large sensi- 
tivities of the overall aerosol simulations to the chosen car- 
bonaceous emission. In this section the results of sensitivity 
experiments are compared to observational datasets, includ- 
ing BC mass measurements and AERONET (Holben et al., 
1998) products. 

Table 4 compares the model with the averaged surface BC 
mass concentrations of the European EMEP (2002-2003) 
and the North American IMPROVE (1995-2001) network. 
The base model underestimates European BC concentrations 
by a factor of 1.6. Decreasing BC emission sizes (SI) im- 
proves the simulation, but surface BC concentrations in Eu- 
rope are still too low by a factor of 1.5. The model also un- 
derestimates surface BC concentrations in North America by 
a factor of 1.8, but with improved correlation coefficients 
compared to the European network. The IMPROVE net- 
work includes mainly rural stations and therefore can better 
be simulated by coarse models, whereas the EMEP network 
includes more urban stations, that observe more detailed sub- 
grid scale features. Comparing surface observations in pol- 
luted regions to a coarse 4°x5° model is questionable, but 
we can conclude that the model underestimates North Amer- 
ican and European surface concentrations. The SI experi- 
ment leads to increased and therefore slightly improved sur- 
face concentrations in those regions. 
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Fig. 3. Changes between BASE and SI (SI -BASE) (upper panels) and (S3 - BASE) (lower panels) for AOT, CDNC [#/cm 3 ], ADE [W/m 2 ] 
and AIE [W/m 2 ]. All maps show differences between the experiments for present day conditions. 


Comparisons between surface concentrations, particularly 
in regions distant from sources, raise the question of how the 
BC is distributed vertically. Figure 2 already demonstrated 
the large variable of BC loads among the single experiments. 


The causes for the large diversity in aerosol load are the dif- 
ferences in aerosol mixing state, solubility and size distri- 
butions, which impacts BC lifetime through wet and dry re- 
moval changes and transport. Recent aircraft measurements 
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Table 4. Comparison between model simulations and EMEP (EU) and IMPROVE (US) measured surface BC concentrations, and AERONET 
absorption optical thickness (AAOT) and optical thickness (AOT). Year 2000 AERONET data are used. Model data are averaged at station 
locations. Gray shadings mark model values that are lower than the observation, red shadings show higher modeled values, and yellow 
shading indicates best simulation. Global and regional comparisons are presented for, North America, South America, Europe, Africa 
(including Arabian peninsula) and Asia. 


EU US AERONET 



BC mass [pg/m 3 ] 



AAOT 





AOT 




EU 

US 

Global 

N. Am. 

S. Am. 

Europe 

Africa 

Asia 

Global 

N. Am. 

S. Am 

Europe 

Africa 

Asia 

# stat. 

14 

55 

77 

33 

9 

13 

10 

6 

77 

33 

9 

13 

10 

6 

Obs. 

0.65 

0.29 

0.018 

0.0061 

0.0014 

0.0035 

0.0032 

0.0031 

0.19 

0.063 

0.011 

0.041 

0.035 

0.031 


BASE 

0.41 

0.16 

0.014 

0.0054 

0.0004 

0.0026 


0.0016 

0.16 

0.057 

0.003 

0.032 


0.017 

SI 

0.44 

0.17 

0.018 

0.0066 

0.0005 

0.0036 

0.0032 

0.0022 

0.19 

0.066 

0.004 

0.040 


0.020 

S3 

0.38 

0.14 

0.010 

0.0032 

0.0003 

0.0018 

0.0032 

0.0013 

0.14 

0.047 

0.004 

0.028 

0.035 

0.015 

S4 

0.31 

0.12 

0.008 

0.0019 

0.0002 

0.0013 


0.0009 

0.11 

0.033 

0.003 

0.022 

0.033 

0.011 


0.0 
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Fig. 4. Global mean ATE and ADE [W/m 2 ] values for all size exper- 
iments, (Sl-4) and the base experiment, BA, for present day condi- 
tions. 


provide BC profile measurements with Single Particle Soot 
absorption Photometers (SP2) onboard NASA and NOAA 
research aircrafts (Schwarz et al, 2006; Slowik et al., 2007). 
In this paper we used the same campaign data and averaging 
technique as in Koch et al. (2009). Figure 5 shows profiles 
measured in North American mid-latitude regions (a, d), in 
the tropics (b,c) and at high latitudes (e-i). The base model 
simulation generally agrees better with some observed pro- 
files in the mid- and high-latitude regions, overestimates BC 
concentrations in the tropics and underestimates BC concen- 
trations in some cases at high latitudes. The sensitivity ex- 
periments show a rather uniform response globally. Smaller 
particle size BC/OC emissions (exp. SI) increase the overall 
BC mass in the troposphere while larger particles (S3, S4) 
leads to a decrease. However, no simulation significantly im- 
proves the vertical distribution of BC mass concentrations, as 
we already get over and underestimations of BC mass con- 
centrations in different latitude regions. Some measurements 


(g and i) are strongly impacted by biomass burning events, 
and these instantaneous events are not reflected in the climate 
model simulation. In summary we dosee some systematic er- 
rors in the vertical distribution of BC (However our regional 
observational coverage is quite limited.): 

- The model overestimates upper-level BC in regions near 
sources, while surface concentrations in those regions 
look reasonable. 

- The model usually, but not always, underestimates sur- 
face BC in the Arctic. 

- The model always fails to reproduce the measured dis- 
continuity (or rapid change) in the vertical distribution 
in the Arctic. 

- Finally, the changes in microphysics do not affect the 
smoothness of the vertical profile, so microphysics does 
not explain the model’s failure to reproduce the ob- 
served sharp discontinuities, however we have to keep 
in mind that we compare observed snap shots to mod- 
eled average mean conditions. 

AERONET (Dubovik et al., 2002) observations provide 
column amounts and absorption of aerosols. These are sea- 
sonally averaged for the year 2000. We use the aerosol op- 
tical thickness (AOT), which is the most reliable measure- 
ment from this network, and the absorption AOT (AAOT), 
the non-scattering part of the aerosol optical thickness, which 
involves a higher degree of uncertainty at low AOT condi- 
tions. 

Figures 6 and 7 show the seasonal comparison between 
the AERONET data and the base model simulations. The 
simulation of AOT shows a good seasonal comparison with 
AERONET. Many seasonal and regional features are well 
captured. Modeled AOT is underestimated in some megac- 
ities, including some European sites and in South America. 
Table 4. Presents the global as well as regional comparisons 
of AOT and AAOT. AOT is best simulated in SI. With larger 
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Fig. 5. Modeled (colored) and observed (black lines) BC mass profiles in [ng/kg] in the tropics and mid-latitudes (a-d) and ARCPAC and 
ARCTAS campaign (e-i). Observations are averaged for the respective campaigns, with standard deviations where available. The two black 
lines in (a) show measurements on two different days. All other profiles give multiple day averages. Mean (solid) and median (dashed) 
observed pro?les are provided except that (g) theARCPAC campaign has distinct profiles for the mean of the 4 flights that probed long-range 
biomass burning plumes (dashed) and mean for the flight that sampled aged Arctic air (solid). 
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Fig. 6. Aerosol optical thickness at 550 nm for clear sky conditions shown for the four seasons. Left column shows model data and the filled 
circles in the right column give the AERONET data. 
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Fig. 7. Aerosol absorption optical thickness at 550 nm for clear sky conditions shown for the four seasons. Left column shows model data 
and in the right column AERONET data. 
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Table 5. Global mean budgets of the BASE experiment, for present day conditions (PD), differences between PD and PI, and differences 
between BASE run and BC reduction experiments under PD conditions. The following variables are listed: Non-biomass burning BC and OC 
emissions, (0 Biomass burning emissions are 3 Tg/a BC and 24.8 Tg/a OC for present-day and 1. T/g BC and 9.1 Tg/a OC for pre-industrial 
times), CDNC (cm 3 ), Liquid water path [mg/m 2 ], total cloud cover [%], the indirect and semi-direct effect (AIE), aerosol direct effect (ADE) 
at the surface and TOA, and total net radiative change (NR) [W/m 2 ]. 



BC1 Emis- 

occ 

CDNC [cm 3 ] 

LWP [mg/m 2 ] 

Cloud Cover [%] 

AIE TO a 

ADE toa 

NRtoa 

ADE surf 


sion 

Emission 





[W/m 2 ] 

[W/m 2 ] 

[W/m 2 ] 

[W/m 2 ] 


[Tg/a] 

[Tg/a] 









BASE (PD) 

4.6 

22.1 

162 

702 


60.1 


-1.78 

-1.75 

-4.62 

BASE 

(PD-PI)® 

0.4 

14.6 

41 

7.0 


0.2 

-0.45 

-0.11 

-0.56 

-1.32 

Mitigation experiments: Delta between experiment and BASE, (EXP-BASE) 

BC 50 (PD) 

-2.3 

0 

4 

-0.2 


-0.05 

0.12 

-0.15 

-0.03 

0.13 

BCFF (PD) 

-3.0 

0 

3 

-0.43 


-0.07 

-0.05 

-0.15 

-0.20 

0.20 

BCOCBF (PD) 

-1.6 

-6.4 

-5 

-2.52 


-0.14 

0.20 

-0.07 

0.13 

0.21 

BC.diesel(PD) 

-1.3 

-0.4 

2 

0.98 


0 

-0.05 

-0.05 

-0.10 

0.11 


BC/OC emission sizes, BASE, S3 and S4 AOT gets further 
underestimated. 

SI shows as well the best simulation of AAOT. Modeled 
AAOT values are systematically too low in South America 
and Asia in all experiments. However, AOT was already too 
low in South America and Asia and therefore the AAOT val- 
ues can’t be evaluated in those two regions. 

The model simulates AOT and AAOT in the strongly dust 
influenced regions of Africa and the Arabian Peninsula best 
in experiments S3. 

In summary, BC surface concentrations are systematically 
underestimated in Europe and the US, but the vertical dis- 
tribution over the North American continent do not show a 
systematic bias. Aerosol optical measures, such as AOT and 
AAOT, are fairly well simulated, and improve with the as- 
sumption of smaller BC particles, which results in larger BC 
loads and more internal mixing. Larger BC particles (exper- 
iments S3 and S4) worsen the simulation. Mineral dust is an 
important absorber and also impacts the comparisonsin this 
study, but we have not carefully examined these effects here. 

5 Black carbon reduction experiments 

Recently, attention has been drawn towards black carbon 
aerosols as a short-term climate warming mitigation can- 
didate. Although the BC direct effect is definitely warm- 
ing, there remains high uncertainty in estimates of indirect 
(cloud) changes associated with reductions, due to the com- 
plex nature of aerosol evolution and its climate interactions. 
Furthermore, BC emission reductions must be achieved by 
controlling individual sources. Therefore co-emitted, cooling 
species would be reduced as well. Here we conduct two ide- 
alized experiments where we simply reduce BC emissions 
from fossil and biofuel burning by 50% (BC_50), and a sce- 
nario where BC fossil fuel emissions are reduced by 100% 


(BC_BCFF). Because OC is not reduced, these experiments 
evaluate the net effect of BC but not of individual sources. 

We also test two more realistic cases where both OC and 
BC from biofuel sources are reduced by 50% (BCOCBF), 
and a scenario where BC and OC emissions from on-road and 
off-road diesel were removed (BC_diesel). Particles emitted 
from diesel engines have the highest BC fraction of any ma- 
jor source, and the OC:BC ratio of 2:1 is common to all en- 
gines. Thus, BC _diesel reflects the likely impact of any con- 
trol strategy that reduces diesel emissions, although it does 
not account for reductions in sulfur emissions which would 
be required to implement advanced controls. On the other 
hand, biofuel sources have different OC:BC ratios, ranging 
from 3:1 or 4:1 for wood cooking stoves to 6:1 for fireplaces 
or cooking with animal waste. The average OC:BC reduc- 
tion of 5.6:1 in the BCOCBF case provides an estimate of 
the response to reductions in high-OC sources, but may not 
exactly reflect the impact of individual sources. See Table 5 
for emission budgets. 

The BC reduction experiments are performed under 
present day conditions and are compared to PD to PI changes 
of the BASE experiment. In addition to the fuel emissions, 
biomass burning emissions for PD runs are 3 Tg/a BC and 
24.8 Tg/a OC. Note that also the non-biomass burning OC 
emissions levels are reduced in experiment BCOCBF and 
BC_diesel. 

The results of all BC reduction experiments are sum- 
marized in Table 5 and Fig. 8. In all experiments reduc- 
ing BC emissions leads to less positive direct aerosol im- 
pacts (ADE), with changes ranging from —0.05 to —0.15 
W/m 2 . BC_50 and BCFF show very similar results for ADE, 
whereas BCOCBF and BC_diesel show much weaker effects. 
This is caused in the diesel case by smaller BC emission 
changes and by a slight decrease in OC emissions. In the 
BCOCBF case some of the negative radiative flux change 
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Fig. 8. Differences in CDNC [cm 3 ], total cloud cover [%], top of the atmosphere cloud (AIE), aerosol (ADE) and surface aerosol radiative 
flux changes [W/m 2 ], between the PD and PI in the BASE case (first panel) and the difference between BASE and the four BC reduction 
cases. Cloud cover and ADE are multiplied by a factor of 10 to fit the color scheme. 


that is achieved by the BC emission reduction is counterbal- 
anced by the positive flux change caused by the reduced co- 
emitted OC emissions. In both cases the reduced OC emis- 
sions over Europe and the US lead to reduced (BC_diesel) or 
reversed (BCOCBF) BC mitigation effects, less positive ra- 
diative flux, in those regions, whereas in Asia and the high 
latitudes both measures would weaken the positive direct 
aerosol flux. 

Table 5 and Fig. 8 also give changes in cloud diagnostics. 
Reducing BC emissions leads to a slight increase in CD- 


NCs in all experiments, except BCOCBF. Changes to CDNC 
come from a change in the aerosol mixing state. Reducing 
the number concentration of BC particles leads to less inter- 
nally mixed BC- sulfate particles, and therefore to an increase 
in the number concentrations of externally mixed sulfate par- 
ticles, which serve very effectively as CDNC. 

First we will discuss the experiments which resulted in 
reduced CDNCs: BC_50, BCFF, and BC_diesel. Region- 
ally, for example over Europe where the strongest change in 
CDNCs occurs, more CDNCs lead to an increase in cloud 
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cover and therefore to a negative cloud effect. However, 
globally, even though there are more CDNCs, there is re- 
duction in cloud cover and positive and negative AIEs rang- 
ing from 0.11 to -0.05 W/m 2 per experiment. These cloud 
cover changes are caused by a competition between indirect 
and semi-direct effects. Fig. 8 shows the surface radiative 
flux changes in the single experiments. Reduced BC loads 
lead to surface heating, due to the missing absorbers in the 
atmosphere. More radiation can reach and warm the Earth 
surface, heat the lower atmosphere and lead to a decrease of 
low-cloud cover and LWP, leading to a positive cloud radia- 
tive effect. These semi-direct cloud changes dominate over 
the smaller cloud changes that result from CDNC enhance- 
ment. 

Case BCOCBF, where BC as well as OC emissions are 
reduced, shows a decrease of -5% in CDNCs. Less CDNCs, 
leads to reduced LWP and cloud cover and therefore logically 
to a smaller (less negative) overall cloud effect, resulting in a 
AIE change of 0.2 W/m 2 compared to the BASE case. 

In Summary: 

- Reducing BC emissions by 50% (BC_50) compared to 
the BASE case leads globally to increased CDNC, de- 
creased cloud cover and LWP and a positive AIE change 
of +0.12 W/m 2 . Together with the direct aerosol impact 
of —0.15 W/m 2 this leads to a net radiative flux differ- 
ence of —0.03 W/m 2 . In principle we would expect a 
negative AIE with increased CDNCs. However we see 
a negative AIE (Fig. 8) in the regions with increased 
CDNCs, e.g. Europe, South East Asia, Easter US, how- 
ever the global mean AIE is positive due to positive AIE 
contributions from remote ocean regions that are not di- 
rectly impacted by CDNC changes. We speculate, that 
those changes, e.g. ADE leads to cooling at the surface 
and TOA over the Southern Hemispheric ocean and pos- 
itive AIE, are caused by semi-direct effects that are re- 
lated to aerosol direct impacts. 

- Removing fossil fuel emissions (BC_FF) leads as well 
to an increase in CDNC, decreased cloud cover but an 
allover negative AIE of —0.05 W/m 2 . Together with the 
direct aerosol impact of —0.15 W/m 2 this leads to net 
radiative flux difference of —0.2 W/m 2 . 

- Removing diesel emissions (BC_diesel) still leads to 
globally increased CDNCs but to a lesser extend than 
the above discussed cases. India and South East Asia 
show a decrease in CDNC. Allover this leads to a strong 
negative AIE over Europe and a global mean AIE of 
—0.05 W/m 2 . Combined with a negative ADE this re- 
sults in a net radiative flux difference of —0.1 W/m 2 . 

- Reducing biofuel sources (BCOCBF) is the only sce- 
nario where we see decreased CDNC concentrations. 
The strong reduction in BC (1.5 Tg/a) and especially 
OC (6.4 Tg/a) particles finally resulted in a reduced 


number concentration of activated particles, which was 
not the case in the prior discussed scenarios, were the 
primary emission reduction (BC and OC) was only 
ranging between 1.7 and 3 Tg/a. In this experiment the 
decrease in primary emission concentration now domi- 
nates over the mixing state effects that lead to localized 
CDNC increases in regions with moderate BC reduc- 
tions. Eventually this leads in the BCOCBF experiment 
to a positive AIE +0.2 W/m 2 and a net radiative change 
flux of +0.13 W/m 2 . 

These model experiments show that the success of BC mit- 
igation in reducing positive radiative fluxesdepends on the 
combination of direct and indirect effects, and the resulting 
change in the net radiative flux. As a reference, our BASE 
model estimates a net radiative flux change of —0.56 W/m 2 
for PD-PI. Our simple BC mitigation scenarios which ignore 
co-emitted species, BC_50 and BCFF, show an decrease in 
the net radiative flux, —0.03 and —0.2 W/m 2 . The two more 
“realistic” scenarios, BCOCBF and BC_diesel, show that 
eliminating diesel emissions will also lead to a decrease in 
the netradiative flux by —0.1 W/m 2 , whereas reducing bio- 
fuel sources would increase the positive aerosol radiative ef- 
fects, due to the combined reduction of BC and OC emissions 
and the resulting reductions in cloud cover. We conclude 
that the amount of the co-emitted species basically decide 
whether BC mitigation strategies are increasing or decreas- 
ing aerosol related radiative fluxes. 

6 Discussion and conclusions 

In this paper we studied the sensitivity of BC effects on ra- 
diative flux changes with a climate model that includes de- 
tailed microphysical processes. Since BA08, MATRIX in- 
ternally and externally mixed aerosol populations are now 
coupled to the radiation scheme and to cloud indirect effects 
(SM10). We examined single processes and particle prop- 
erties and their impacts on cloud and aerosol radiative flux 
changes. This study allows us to evaluatewhich processes 
need to be further studied, and what observational data are 
needed. 

First we discussed the importance of initial size distribu- 
tions for particulate emissions such as BC and OC. Evaluat- 
ing the model with ambient BC mass and AERONET AOT 
and AAOT measurements leads to the conclusion that best 
results are seen when BC/OC emission sizes are assumed to 
be between 0.01 and 0.05 pm. Aerosol microphysical mod- 
els such as MATRIX now require particulate emission inven- 
tories that include information about aerosol mass, number, 
surface area, composition, and mixing state, and possibly in- 
cluding subgrid scale effects, such as source-plume mixing, 
that can not be treated by current climate models. 

Koch et al. [2009] showed that the previous generation 
of aerosol models without aerosol microphysics and mixing 
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state influencing optical properties generally underestimate 
AAOT. In contrast, this study demonstrated a greatly im- 
proved agreement with retrieved AAOT, as the internal mix- 
ing of BC enhances absorption. On the other hand, Koch 
et al. (2009) found that these older models generally did not 
underestimate BC surface concentration; the rapid BC ag- 
ing of the microphysical scheme in our model caused faster 
washout near source regions and consequent underestimation 
of BC surface concentrations. Thus it seems that the micro- 
physics required to improve BC optical properties has ex- 
cessively reduced model BC surface concentrations. Careful 
examination of regional tendencies may help resolve the new 
discrepancy. 

However in order to better understand this behavior we 
will need to better evaluate the aerosol mixing state simulated 
by our model. The largest volume fraction is taken up by 
aerosol water, a quantity that is hardly validated in aerosol 
models. Chemical and morphological analysis of different 
individual particles is needed, along with size selective bulk 
analysis, to understand such processes, taking into account 
aerosol microphysical processes. 

One possible approach to observe mixing state includ- 
ing BC coating estimates was discussed by Shiraiwa et 
al. (2008), where the mixing state of BC in Asian air 
masses was examined by a combination of aerosol mass 
spectroscopy (AMS) and single-particle soot photometer SP2 
data. The AMS data provide a detailed chemical analysis 
of the measured species, and the SP2 was able to measure 
the size distribution, and the mixing state of BC. Shiraiwa 
et al. (2008) found a median value of the shell/core diame- 
ter ratio increased to 1.6 in Asian continental and maritime 
air masses with a core diameter of 200 nm, while in free tro- 
pospheric and in Japanese air masses it was 1.3-1. 4. It was 
estimated that internal mixing enhanced the BC absorption 
by a factor of 1.5-1. 6 compared to external mixing. The cal- 
culated absorption coefficient was 2-3 times higher in Asian 
continental air masses than in clean air. A further approach 
would be the combination of AMS with single particle mass 
spectroscopy ATOF-MS (Prather et al, 1994), or chemi- 
cal ionization mass spectroscopy, CIMS (Hearn and Smith, 
2004) data. Our ongoing research will focus on the valida- 
tion of mixing state with these newly emerging datasets. 

In light of the analyzed model sensitivities towards mi- 
crophysical BC characterization, we performed a set of BC 
mitigation studies. The model suggests that BC mitigation 
is always beneficial if only BC sources are reduced, as the 
direct radiative flux is lowered, and the effects on clouds are 
weak. BC mitigation led to enhanced CDNC production and 
a slightly stronger indirect effect; this finding is in contrast to 
assumptions about reducing primary particle concentrations. 
However, semi-direct effects, induced by changed vertical 
heating profiles in the atmosphere, as well as the first and sec- 
ond indirect effect, are very important components of there- 
sponse to altered BC emissions. The two more ‘realistic’ 
scenarios which include reduction of OC show that mitigat- 


ing diesel emissions will decrease aerosol radiative fluxes by 
-0.1W/m 2 , whereas reducing biofuel sources could even in- 
crease radiative fluxes, mainly due to the impacts on clouds. 
We can not quantify the forcing uncertainties of the BC re- 
duction experiments, but we want to point out that in light of 
the large uncertainties already involved with the direct effect, 
as we discussed in this paper, and the vast complexities of the 
indirect aerosol effects the uncertainty of the here presented 
numbers is large. 

This study did not include ice cloud feedbacks or BC- 
snow ice albedo feedbacks (Koch et al. 2009) that may en- 
hance warming by BC. Despite the various limitations ex- 
plored in this study regarding emission sizes and mixing state 
of BC particles, our model results suggest that a reduction 
in BC diesel emissions would help reduce positive forcing 
even when including cloud changes. However, we stress that 
strategies that include a reduction in BC emissions should 
not delay any GHG reduction plans, as the likely reductions 
in positive radiative fluxes from BC emissions are small and 
uncertain (from +0.12 up to -0.19W/m 2 ) compared to GHG 
forcings of 2.7 W/m 2 and growing. Furthermore, BC reduc- 
tions also reduce co-emitted species, many of which have 
negative forcings, such as OC and sulfate. For sources suf- 
ficiently richin black carbon, mitigation is beneficial, but 
aerosol microphysical processes and their interactions with 
atmospheric dynamics need to be better constrained in order 
to estimate climate warming benefits from those measures. 
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